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ABSTRACT 


In  Section  I it  is  shown  by  using  the  simplest  reflecting 
surf'ace--an  infinite,  pressure  release  plane  surface-- 
that  at  least  the  Fresnel  phase  approximation  and  a reasonable 
directivity  pattern  approximation  are  necessary  before  accurate 
predictions  in  surface  reflection  and  scattering  problems  can 
be  made.  The  average  Rayleigh  reflection  coefficient  for  random 
penetrable  surfaces  of  various  roughnesses  is  calculated  in 
Section  II.  In  Section  III  an  integral  expression  is  derived 
for  the  acoustic  field  of  a point  source  in  a surface  duct  with 
a rough  boundary.  The  integral  was  obtained  by  Green's  functions 
and  includes  the  rough  surface  condition  by  means  of  an  effective 
impedance  condition  which  also  depends  on  the  properties  of  the 
medium. 
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I.  THE  REFLECTION  OF  AN  ACOUSTIC  BEAM  BY  AN  INFINITE, 
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ATTENUATING  BOTTOMS 

III.  ACOUSTIC  FIELD  IN  A SURFACE  DUCT  WITH  A ROUGH  BOUNDARY 

A.  Introduction 

B.  Green's  Function  Solution 


REFERENCES 


v 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


I 


I.  THE  REFLECTION  OF  AN  ACOUSTIC  BEAM 
BY  AN  INFINITE,  PRESSURE  RELEASE  PLANE  SURFACE 

The  prediction  of  the  field  scattered  by  a known  surface 
configuration  or,  conversely,  the  prediction  of  the  surface  configu- 
ration from  its  scattered  field  is  a difficult  and  complex  problem. 

Due  to  the  extreme  complexity  of  the  problem,  all  of  the  solutions  of 
the  integral  equations  obtained  to  date  have  contained  numerous  approxi- 
mations. One  of  the  fundamental  difficulties  associated  with  surface 
scattering  has  been  the  determination  of  the  ranges  of  validity  and  the 
interaction  of  the  various  approximations. 

I I 

The  usual  approximations  can  be  put  into  five  categories: 

(l)  formulation,  (2)  boundary  values,  (3)  source  directivity, 

(4)  phase  approximation,  and  (5)  surface  representation.  It  was 
shown  in  the  Final  Report  under  Contract  N00024-69-C-1275  (April  1970) 
that  the  potential  formulation  of  reflection  and  scattering  problems 
has  distinct  advantages  over  the  Helmholtz'  and  Green's  functions 
formulations  and,  in  the  limit  of  a plane  surface,  reduces  to  the 
Rayleigh -Sommer f eld  formula.  The  boundary  value  for  a plane,  pressure 
release  surface  is  known  exactly.  The  boundary  values  for  arbitrary 
penetrable  surfaces  will  be  the  subject  of  later  progress  reports. 

i 

Surface  representations  for  randomly  rough  surfaces  will  also  be 
treated  further  in  future  reports. 

I 

' 

At  this  time  only  an  infinite,  pressure  release  plane  surface 
will  be  considered,  since  this  is  the  simplest  of  surfaces  and  the 
results  are  easily  verified  by  the  method  of  images.  By  treating  a 
plane  surface,  attention  can  be  focused  exclusively  on  the  phase  and 


source  directivity  approximations  that  will  be  necessary  in  the 
treatment  of  arbitrary  penetrable  surfaces  and  the  inverse  problem, 
and  certain  misconceptions  about  the  phase  and  directivity  approxima- 
tions can  be  exposed. 

To  gain  insight  into  what  approximations  might  be  necessary,  it 
is  useful  to  recall  two  experimentally  established  properties  of  the 
intensity  reflected  by  a pressure  release  plane  surface:  (l)  the 

range  dependence  is  that  predicted  by  geometrical  acoustics,  namely 
l/(r  +r^),  and  (2)  the  results  predicted  by  the  image  technique 
are  indeed  obtained.  Drawings  AS-69-^78  and  AS-69-^79  demonstrate 
these  properties.  In  light  of  the  experimental  facts,  two  conclusions 
are  inescapable:  (l)  the  Fraunhofer  phase  approximation  is  never 

correct  when  one  is  concerned  with  surface  reflection  or  scattering 
since  it  predicts  a 1/r^r.^  spreading  loss;  and  (2)  the  analytic 
solution  should  reduce  to  the  image  solution.  The  first  conclusion 
applies  to  arbitrarily  rough  surfaces  since  in  the  limit  as  the  rough- 
ness approaches  zero,  the  results  must  reduce  to  a plane  surface. 

Hence,  second  and  possibly  higher  order  terms  must  be  used  in  the 
expansion  of  r^  and  r^  in  the  phase  (Fresnel  approximation);  otherwise, 
the  correct  range  dependence  will  never  be  obtained. 

In  previous  reports  (see,  for  example,  Final  Report  under 
Contract  N00024-70-C-1275) > it  has  been  shown  that  the  pressure 
reflected  from  a plane  pressure  release  surface  to  some  receiver  at 
A is  given  by: 
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where 


1)  if  the  amplitude  factor  is  suppressed,  Pq  is  simply  the 
directivity  pattern  of  the  source, 

2)  r and  r,  are  the  distances  from  the  source  and  the  receiver 

o 1 

respectively  to  the  surface  element  dS, 

3)  n is  the  unit  normal  out  of  the  surface  at  the  point  enclosed 
by  dS, 

4)  is  the  unit  vector  pointing  from  the  surface  element  dS 
along  r^  to  the  receiver,  and 

5)  S'  represents  the  insonified  surface  area. 

The  only  assumption  contained  in  this  result  is  kr»l.  The  geometry 
is  depicted  in  the  accompanying  diagram. 


The  origin  is  fixed  at  the  point  where  the  axis  of  the  acoustic  beam 
intersects  the  plane.  From  simple  trigonometry,  it  can  be  shown  for 
a plane: 
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In  general,  Eq.  (l)  cannot  be  evaluated  without  making  certain 
approximations.  Previously,  it  was  pointed  out  that  at  least  the 
Fresnel  approximation  was  required  in  the  phase.  The  Fresnel  phase 
approximation  is  given  by  the  equation 

1 2 12 

ik(r  +r  ) = ik(rn  + r + ax  + — x^  + — y“)  + higher  order  terms  , 
1 o lo  oo  R R D ’ 

1 (4) 
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Elsewhere  r and  r.,  will  be  replaced  by  r and  r,  . 
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All  that  remains  now  is  to  find  a suitable  approximation  to  the 
directivity  pattern.  For  a nonisotropic  source  there  are  essentially 
two  tractable  approximations  for  the  directivity  pattern.  The  first 
is  the  simple  aperture  type  function  where  it  is  assumed  the  insoni- 
fication  is  a constant  over  the  region  of  the  surface  defined  by  the 

-3  cLB  points  of  the  source  and  zero  elsewhere  on  the  surface.  The 
geometry  is  shown  in  Dwg.  AS-69-635.  If  the  beamwidth  of  the  source 

is  $ deg,  then  the  following  expressions  hold: 


= r^sin(  0/2)/sin(  0^+<J>/2)  , 

(5) 

A = r sin(®/'2)/sin( 9.  -4>/2)  , 

(6) 

l = r tan( $/2)  , and 

(7) 
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Then : 


K = 3./20  log  e 
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where 

1)  K is  a constant  which  provides  the  -3  dB  value  at  the 
appropriate  point  off  the  beam  axis, 

2)  a is  the  semimajor  axis  of  the  insonified  ellipse,  and 

3)  P is  the  semiminor  axis. 

Note  that  the  effect  of  the  range  differences  from  one  point  on  the 
ellipse  to  another  has  been  ignored  in  setting  the  -3  dB  points.  These 
expressions  are  calculated  in  detail  in  the  Quarterly  Progress  Keport 
’under  Contract  N00024-69-C -1275  for  April  - June  1969.  In  the  calcula- 
tions, a beamwidth,  $>,  of  8 deg  will  be  assumed  since  this  is  the 
beamwidth  of  the  source  used  in  the  experimental  program. 


If  now  the  ellipse  is  approximated  by  a rectangle  of  dimensions 
2 i,  A.+A^,  then  the  pressure  reflected  by  a plane  surface  where  the 
directivity  pattern  is  assumed  to  be  a simple  aperture  may  be  calculated. 
Using  Eqs.  (2),  (3),  and  (4)  in  Eq.  (l)  gives 
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The  integrals  in  Eq.  (12)  are  by  definition  Fresnel  integrals  and  upon 
forming  the  product  of  p(A)  and  its  complex  conjugate,  the  intensity  is 
obtained: 


Ip  = p*(A)p(A) 


where 


and 


dM 


(15) 


Then 
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It  is,  however,  possible  to  make  a considerably  better  approximation 
to  the  directivity  pattern.  If  the  side  lobes  may  be  ignored,  a 
suitable  functional  form  is 


P 

o 


exp[-Kty2/(t>/2)^l 


(22) 


where 

1)  K and  $ have  been  previously  defined,  and 

2)  y is  the  angular  displacement  measured  from  the  beam  axis. 

Since  the  integration  is  to  be  performed  over  the  surface,  it  is 
necessary  to  express  Eq.  (22)  as  a function  of  the  surface  variables  x 
and  y.  A useful,  analytically  tractable  approximation  to  Eq.  (22)  is 


P = 
o 


exp 


(£3) 


2 3 

This  lorm  was  suggested  by  Clay,  and  Horton  and  Muir.  The  directivity 
pattern  given  by  Eq.  (23)  represents  a single  lobe*  source  whose  surface 
insonification  is  an  ellipse  of  semimajor  axis  ot,  and  semiminor  axis  p. 
With  the  introduction  of  Eq.  (23)  the  integration  limits  may  now  be 
taken  as  infinite  for  most  practical  cases.  This  consideration  can  be 
very  important  in  connection  with  arbitrary  surfaces,  since  often  the 
integrability  of  the  equations  will  depend  on  the  extent  of  the  limits. 


It  is  possible  to  generalize  Eq.  (23)  to  a multi-lobe  directivity 
pattern;  for  an  N-lobe  source  the  equation  becomes 


i - 1 


10 


In  Dwg.  AS-7O-8OO,  the  experimentally  measured  directivity 
pattern  of  the  100  kHz  source  used  in  scattering  experiments  is  com- 
pared with  Eq.  (22).  Drawing  AS-7O-809  compares  Eqs.  (22)  and  (2j) 
for  various  source  grazing  angles 

Now  the  intensity  reflected  from  a plane,  pressure  release  surface 
can  he  calculated.  Substituting  Eqs.  (b)  and  (23)  into  Eq.  (l)  and 
taking  the  product  of  p(A)  with  its  complex  conjugate  gives 


—00 


(2b) 


ka( x-x ' ) 


dxdydx ' dy ' 


Make  the  transformations  5 = x-x',  7 = y-y' , x = (x+x')  2,  and 
y = (y+y')  2 then  the  integrations  follow  easily,  yielding 


12  2 ” 

exp  -KRi,  a"/ 2af  ( 1+p) 

T„  ^ 2 — , » ( 0 

(roo+rlo)(roo+  rio  Sin'  V Sin^r}  V(  1+PX  1+Q) 

where 

P - (kri  ka  ) , and 
Q (lCR/kp  ) . 

f and  Q are  approximately  zero  for  the  case  of  interest  here.  (This  is 
merely  a restatement  of  kr»]  . ) 
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COMPARISON  OF  MEASURED  DIRECTIVITY  PATTERN 
AND  ANALYTIC  DIRECTIVITY  PATTERN 
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The  solution  obtained  using  image  techniques  is 


j image 
P 


( ro+rl )' 


exp[-2Ki|<2/(i)/2)2. 


(ro+ri)' 


(26) 


where  i is  now  the  angle  to 

the  receiver  measured  from  the 

beam  axis  of  the  image  source . 

The  geometry  is  given  in  the 

accompanying  diagram.  Through 

simple  geometry,  y is  found  in 

terms  of  r , r,  , 0.,  and  9 : 
oo  lo’  3.  r 


\jr  = Arctan 


rn  sin(  0. ~9  ) 
lo  i r 

r +rl  cos ( 0. -6  ) 
oo  lo  i r _ 


(27) 


and  by  the  Law  of  Cosines 


r +r.  +2r  r,  cos (0.-0  ) 
oo  lo  oo  lo  i r 


(28) 


In  the  specular  direction,  examination  of  Eqs.  (25 ) and  (26)  reveals 
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T image 
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(r  +r.  V 

oo  lo 


(29) 


The  image  solution  Eq.  (26)  is  compared  with  Eqs.  (13)  and  (25) 
for  various  source  grazing  angles  and  receiver  ranges  in  Dwgs.  AS-70-801, 
AS-70-802,  and  AS -70 -803.  The  Fresnel  solution  with  the  realistic 
directivity  pattern  (Modified  Fresnel  solution)  is  clearly  superior  to 
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the  Fresnel  solution  with  an  aperture  type  directivity  pattern  (Fresnel 

solution) . It  is  interesting  to  note  that  at  low  grazing  angles  and 

for  receiver  ranges  of  the  same  magnitude  as  the  source  range,  the 

distortion  in  the  Modified  Fresnel  solution  corresponds  closely  to  the 

distortion  in  the  directivity  pattern  (Eq.  23),  hut  is  reflected.  As 

r,  /r  increases,  the  distortion  in  the  Modified  Fresnel  solution 
lo'  00 

decreases  until  for  large  r.  /r  the  Modified  Fresnel  is  nearly  equal 

lo'  00 

to  the  image  solution. 
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II.  CALCULATION  OF  REFLECTION  COEFFICIENTS 
FOR  ROUGH  ATTENUATING  BOTTOMS 

In  this  section  the  Rayleigh  reflection  coefficient  is  modified 

4 

to  account  for  bottom  roughness.  The  method  of  Mackenzie  is  followed 
so  that  the  effects  of  roughness  will  be  easily  observable.  From 
Dwg.  AS -70 -795  it  is  obvious  that  to  include  the  effects  of  roughness 
it  is  only  necessary  to  express  the  local  incident  and  refracted  angles, 
$ and  in  terms  of  the  angles  from  the  mean  plane,  0,  and  0O,  at 

each  point  on  the  surface.  These  are  related  through  the  angle  p which, 
for  a randomly  rough  surface,  will  be  a stochastic  variable.  From 
the  drawing  it  is  clear  that 

+ P , 0=0+3  . (50) 

It  is  assumed  that  at  any  point  on  the  surface,  the  reflection 
coefficient  is  given  by  the  usual  Rayleigh  coefficient  referenced  to 

4 

the  tangent  plane.  Following  Mackenzie,  this  expression  is 


’•.Vie  re 
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(31) 


P1c1sinJ'2 

P2c2sin01 


p c ,sin0  - p c,  rinO 

1 l : 

p. c^sinO^  + p^c^sinO 


R 


(33) 


Using  Eq.  (30)  in  Eq.  (33)  and  the  appropriate  trigonometric  identities 
gives 


pQc,.  ( sinQ^cosP  + cos0^sinf3)  - p^c^(  sin©„cosp  + cosfiLsinp) 
p.  c ( sinfl  cosp  + cos0^sinp)  + ( sin0,.cosp  + cosS^sinP) 


q = tanP 


so  that  q is  the  slope  at  the  given  point  on  the  surface,  and  let 


A_  = poc2sin01  - p^c^sin©^ 


A+  = p^c.^sinQ^  + p1c1sin0.. 


B = p0c0cos9,  - p.c,cos9_  , and 
1 1 1 d 


B+  = 020^0030^  + PqC^COS0? 


This  gives 


A + qB  (A  + qB  ) 


A+  + 


1 + — ^ 
^ + 


Expanding  this  in  terms  of  q gives 


A /A  B - A B 


/ A"1  K 
|(-1  \V  v 


where  A_/A+  is  the  Rayleigh  coefficient  for  a pl-uie  surface. 


The  parameter  S may  now  he  redefined  in  terms  of  the  angles  6 and 

6 as 


p c sin©2 
o2c2sinei 


After  a little  algebraic  manipulation  the  Rayleigh  coefficient  becomes 


l!!Vf£i>y'(-i) 

tan0^  [i+s '] 


K K I,  „ 

q <tan0 


1+r — a~  s' 

tan6? 


• (to) 


Since  q is  a random  variable,  to  find  <R>  it  is  only  necessary  to 
replace  q by  <q>  in  the  series  above.  It  is  assumed  that  t}  has  a zero 

ciK+1. 

mean  Gaussian  distribution.  This  implies  that  all  odd  moments  <r\  > 

are  zero. 


To  include  the  effects  of  attenuation  in  the  bottom  it  is  now 


necessary  to  allow  S'  to  be  complex: 


S'  = x + iy 


The  quantities  x and  y are  related  to  Mackenzie's  parameters  g and  h by 


. . h+ig 
x + iy  = 


where  o is  the  impedance  ratio  at  the  boundary.  If  only  the  first 
nonzero  term  of  the  series  is  retained,  <R>  may  be  put  in  the  form 


^ “ C'+iD' 
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where 


A'  = 


[(1+x)  -y  l(l-x)  + :y  (1+x)  + 2<r,2>f!!fi  _±) 
L J ytane  J 


tanei  , 2 2s 

x + u^e  - (x  -y  } 


B’  = 


C'  = 


D’  = 


r 1 /tan 0 

-y[(l+x)^.y-]+  2y  (1-x-)  + ^ 


tan©  j 

y + i^r  xyl 


[(l+x)  -yj(  1+x)  - 2y"(l+x) 
y £(  1+x)'  -y  J + 2y(  1+x) 


. and 


(kh) 


Finally,  the  magnitude  of  the  reflection  coefficient  including 
attenuation  and  roughness  is  given  By 

2 

| <K>|  = ( — rr ) |(A'C'+B'D’)  + (b'C'-A'D')  J . (^5) 

'C'+D'/  L J 

In  calculating  the  RMS  .-lope  <t  >,  it  was  assumed  that  the  surface  was 
shadowed  at  all  points  where  the  slope  was  less  than  the  tangent  of 
the  incident  angle.  This  assumption  results  in  a slmpl<  :hang<  Ln  th< 
limits  of  the  integral  that  defines  <1^  and  gives  a value  of  > that 
significantly  differs  from  the  true  RMS  slope  only  at  small  grazing 
angles . 

Three  separate  fluid  bottoms  with  different  densities  and 
velocities  were  considered.  These  correspond  to  three  of  the  samples 
used  by  Mackenzie.  The  calculated  quantity  is  the  magnitude  of  the 
reflects  i :oefficien1  a<  given  by  Eq.  This  equate  n toes  not 


an  equation  for  the  full  scattering  coefficient.  It  simply  represents 
the  change  in  the  reflection  coefficient  resulting  from  the  introduction 
of  random  roughness  into  the  geometry  of  the  problem . In  subsequent 
reports  the  angular  distribution  will  be  combined  with  these  results  to 
give  an  equation  for  the  scattering  coefficient  for  penetrable  rough 


surfaces . 


All  of  the  calculated  curves  were  evaluated  at  5-0  kHz  for  ease  of 
comparison  with  Mackenzie's  data.  Several  different  values  of  RMS  slope 
were  calculated,  but  only  the  one  which  produced  the  largest  variations 
in  the  reflection  coefficient  was  plotted.  Drawings  AS-70-796  and 
AS-70-797  show  cases  where  c^/c  <1.  For  these  cases  it  is  possible  to 
have  an  angle  of  intromission;  however,  neither  of  the  curves  shown  exhibit 
this  behavior  because  of  the  assumed  attenuation.  Drawing  AS-70-798 
shows  a case  where  c.,  c >1.  Again,  the  critical  angle  behavior  of  this 
bottom  is  absent  due  to  the  attenuation.  It  should  be  noted  that  where 
c0/c  <1.  the  reflection  coefficient  that  includes  roughness  lies,  for 
the  most  part,  below  the  plane  reflection  coefficient,  while  for 
Cg/c  >1,  the  opposite  is  true. 

In  conclusion,  the  overall  effect  of  including  roughness  in  the 
reflection  coefficient  is  very  slight.  In  none  of  the  cases  considered 
was  the  variation  more  than  about  2 dB  for  the  roughest  bottoms  and. 
indeed  in  most  cases,  the  variation  was  only  a few  tenths  of  a decibel. 
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REFLECTION  COEFFICIENT  AS  A FUNCTION  OF  GRAZING 
ANGLE  FOR  A CLAYEY  FINE  SILT  BOTTOM 
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REFLECTION  COEFFICIENT  AS  A FUNCTION  OF  GRAZING 
ANGLE  FOR  A SAND  - SILT  - CLAY  BOTTOM 


SILTY  VERY  FINE  SAND  - 5.0  kHz 
MACKENZIE  SAMPLE  No.  28 
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REFLECTION  COEFFICIENT  AS  A FUNCTION  OF  GRAZING 
ANGLE  FOR  A SILTY  VERY  FINE  SAND  BOTTOM 


III.  ACOUSTIC  FIELD  IN  A SURFACE  DUCT 
WITH  A ROUGH  BOUNDARY 

A.  Introduction 


The  propagation  of  acoustic  waves  in  a surface  duct  with  an 

irregular  "boundary  is  treated  by  means  of  Green's  functions.  Ray 

theory  has  been  applied  to  this  problem  by  various  authors / ' but 

for  shallow  water,  low-frequency  propagation  ray  theory  is  inadequate. 

Normal  mode  wave  theory  has  also  been  used  ^ in  examining  this 

problem.  A useful  survey  of  both  ray  theory  and  wave  approaches  has 

23 

recently  been  given  by  Schulkm. 


The  following  treatment  will  be  very  similar  to  that  of  Bucker's 

but  will  use  Green's  functions  and  a different  impedance  condition  in 

the  basic  integral.  The  impedance  condition  now  depends  ori  the  nature 

of  the  medium  as  well  as  the  random  boundary.  A rough  boundary  can  be 

approximated  by  a flat  boundary  with  some  effective  impedance  (or 

reflection  coefficient)  in  the  case  of  wave  guide  propagation,  but 

the  expression  for  this  impedance  (or  reflection  coefficient)  will  be 

completely  different  from  the  forms  for  reflection  from  a single 
14,17,19,20 

uneven  surlace. 


When  the  rough  surface  is  replaced  by  an  impedance  condition 

(or  effective  reflection  coefficient),  the  propagation  problem  can  be 

22 

solved  by  classical  methods.  As  Bucker  indicates,  this  approach  is 
useful  for  the  case  where  both  source  and  receiver  are  in  the  duct. 
This  case  is  referred  to  as  first-order  scattering  (scattering  out  of 
a mode).  When  the  source  is  in  the  duct  and  the  receiver  is  below  the 
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duct,  then  one  must  consider  second-order  scattering  (scattering  into 

a mode).  The  treatments  of  second-order  scattering  given  by  Schweitzer"^ 

8 

and  Van  Ness  are  believed  to  be  in  error.  The  present  treatment  will 
only  consider  the  first-order  effects  of  the  rough  boundary. 

B.  Green's  Function  Solution 


The  wave  equation  for  a point  monopole  source  of  unit  strength, 
with  angular  frequency  cn,  is  given  by 


p 

A p 


-4—  A = 5 (?-?  )e'iut 

c^(z)  3t 


(^6) 


where 

1)  p represents  the  pressure, 

2)  c(z)  is  the  sound  velocity  (variable  in  the  coordinate  z),  and 

3)  6(r-r  ) is  the  three  dimensional  Dirac  delta  function. 

If  cylindrical  coordinates  (r,  9,  z)  with  azimuthal  symmetry  are 
assumed  and  the  time  factor  of  exp(-ia't)  is  suppressed,  then 


- — 6(  r ) 5( z-z  ) 
r o 


where  l)  k(z)  = cn/c(z)  represents  the  wave  number,  and  2)  the  source 
is  located  at  z=z^  and  r=o,  as  indicated  in  Dwg.  AS-7O-767.  The 
vertical  depth  coordinate  z varies  from  o S z g «,  and  the  range 
coordinate  r varies  from  o < r < °°. 
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The  pressure  is  then  exp(-int)  times  the  Green's  function  G of 
Eq.  (47).  The  "boundary  conditions  imposed  on  Eq.  (kj)  are 

a)  G must  satisfy  a radiation  condition  for  r — * 00  and  z — * », 


"b)  G 


, dG 


= 7 


m 


where  y represents  the  impedance  (or  corresponding  effective  reflection 
coefficient)  condition  found  "by  the  prescription  given  by  Lysanov."0 


In  the  Green's  function  approach,  Eq.  (47)  is  separated  into  the 
following  forms: 


[Lr+\ir]  GL(r,\i)  = -26(r)  , (1+9) 

CLz-\]  G2(z,zo,-\1)  = -6(z-zq)  , (50) 

where  the  differential  operators  are  defined  by 

Lr  = (r  fe)  ’ and  (51) 

d2  2 

Lz  = — J + k^(z)  > and  (52) 

dz 


where  is  the  separation  constant.  G and  G,,  are  defined  as  the 
resolvent  Green's  functions  when  they  satisfy  Eq.  (49)  and  Eq.  (50), 
respectively,  and  are  subject  to  the  previously  stated  boundary  condi- 
tions. The  separation  constant  is  a complex  parameter  so  restricted  as 
to  assure  a unique  solution  for  Eqs.  (49-50). 
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The  Green's  function  G of  Eq.  (kj)  is  now  given  by  a complex 
convolution  of  G and  G,  : 

0(r,z,zQ)  = — Jc  G1(rA1)  G2(z,zo,-\1)  dXL  , {lyj) 

where  the  contour  C separates  the  singularities  of  G,  and  G„  and  can 

I CL 

be  closed  around  those  singularities  of  either  one  which  allow  a 
uniqueness  condition  to  be  satisfied  for  G and  G,- . 

The  solutions  to  Eq.  (U9)  and  Eq.  (50)  that  satisfy  the  required 
boundary  conditions  will  be  given  now.  For  the  exp(-invt)  time  factor, 
these  solutions  are  given  by 


G1(r,\i)  = inH^  (|r)(|  = \J'Oo  < arg  ^ < 2n ) 


and 


G„( z , z ,-\) 
c-  0 1 


n (z,|)[n  (z  ,i)  - X-n  (z  , |) 
1 *-  _ o 10 


W(n. ,n^ ) 


z < Z < 00 

o 


. (55) 


n (z  , |)[n,.(z»t)  - X”n  (r  , 01 

G ,(  Z , Z ,-X.)  = — rri j ; o < Z < Z , ( r>6) 

■ o 1 W(n  ,ru  ) o 


where  H^( |r)  is  the  Hankel  function  of  the  first  kind,  and 


n .( °>  0 - in'  (o,0 

ti,(o,0  ~ xn'  (o,0 


( 57) 


The  Wronskian  of  n and  n ^ is  represented  by  W(n  .n^).  where  is  a 
solution  to  the  homogeneous  part  of  K.j . ( ‘ >0 ) which  has  outgoing  waves 


at  z -> +■»,  and  which  has  outgoing  waves  at  z -»  -°°.  The  primes 
indicate  differentiation  with  respect  to  z. 

When  Eq,  (5*0  and  Eq.  (55)  are  used  in  Eq.  (53),  the  Green's 
function  for  Eq.  (*+7)  is  obtained: 


G(r,z.zQ)  = 


L 


il(z<S)[n2(z>,!)  -X~n1(z>, ?)]  H^(£r)|dg 
W(n2,n1) 


(58) 


where  z^  and  z>  denote  the  smaller  or  larger,  respectively,  of  the 

variables  z and  z . 

o 

When  Eq.  (58)  is  integrated  by  Cauchy's  residue  theorem,  the  normal 
modes  (plus  any  branch  line  integrals)  are  obtained.  The  poles  yielding 
the  normal  modes  are  given  by  the  zeroes  of  the  expression 

jn^o,  |)  - rn|  (o,g)J  • 

Equation  (58)  represents  the  desired  formula  for  the  propagation 
of  sound  in  a surface  channel  with  a rough  boundary.  It  should  be 
noted,  however,  that  because  of  the  formulation,  only  first-order 
effects  of  the  rough  surface  are  considered. 

Future  work  will  deal  with  the  application  of  Eq.  (58)  to 

specific  velocity-depth  profiles,  the  first  being  the  bilinear  profile. 

22 

The  results  will  then  be  compared  with  Bucker's,  who  also  used  the 
bilinear  profile  with  a different  interpretation  for  the  impedance 
condition . 
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